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Abstract. A field theoretical renormalization group approach at two loop level is applied to the homoge- 
neous spin-1 Bose gas in order to investigate the order of the phase transition. The beta function of the 
system with d — 4 — e dimensions is determined up to the third power of the coupling constants and the 
system's free energy on the border of the classical stability is given in next to leading order. It is found 
that the phase transition of the interacting spin-1 Bose gases with weak spin-dependent coupling constant 
values is of first order. 

PACS. 03.75.Mn Multicomponent condensates; spinor condensates - 05.10.Cc Renormalization group 
methods - 05.70.Fh Phase transitions: general studies - 64.60.-i General studies of phase transitions 



1 Introduction 

Bosc-Einstcin condensation (BEC) of dilute, interacting, 
scalar atomic gases is widely believed to be a continuous 
phase transition. Mean-field theory results based on per- 
turbation theory are however contradicting. The simplest 
approximations call for a continuous phase transition [1] , 
while according to the more sophisticated Hartree-Fock 
(Popov) approximation BEC is of first order [2-4]. The 
contradicting results are understood to be the consequence 
of critical phenomena. Namely, close to the critical point, 
the different Feynman diagrams develop to the same mag- 
nitude, and a perturbative treatment becomes inaccurate. 
The different renormalization group calculations are indi- 
cating that BEC is a continuous phase transition falling 
into the same universality class as the 0(2) model of field 
■ theory [5-8]. 

In experiments made with dilute gases of alkali atoms 
in optical traps [9-14] the particles have internal spin de- 
grees of freedom. Such systems at low temperatures can 
| be modeled by Hamiltonians with multiple coupling con- 
, stant interactions [15-17] in the s-wave scattering limit. 
For such a situation, with multiple coupling constants, 
(and in the homogeneous case) the order of the transition 
is not necessarily of continuous type, as e.g. in the case of 
the field theoretical <f> model with cubic anisotropy, where 
fluctuations can induce the transition to be of first order 
[18,19]. _ 

In this paper we treat the problem of the homogeneous, 
spin-1 Bose gas, where two coupling constants arise natu- 
rally in the low energy limit. In the absence of an external 
magnetic field depending on the magnitude of these two 
parameters two possible Bose-Einstein condensed phases 
exist, namely, the ferromagnetic phase, in which the sys- 
tem favors a macroscopic magnetization and the polar 
phase, in which the system prefers no magnetization. The 



gas condensing to the ferromagnetic phase is called the 
ferromagnetic gas, while the other is the polar gas (see 
Eqs. (4) and the discussion below). We assume zero ex- 
ternal magnetic field. For such a system mean-field theory 
results are also contradicting. In the Hartree-Fock approx- 
imation [20] the phase transition of both the polar and 
ferromagnetic gases is of first order. The jump in parame- 
ters is the function of the bigger coupling constant, which 
is responsible for non spin-flip scatterings. However this 
strongly first order type of transition is considered to be 
an artifact of the Hartree-Fock approximation such as in 
the case of the scalar Bose gas. On the other hand the 
Hartree approximation [21], which is in a way a simpler 
mean field approximation yields a continuous Bose con- 
densation in the polar case, while a first order one in the 
case of a ferromagnetic Bose gas. In the latter case the 
jump is a function of the smaller coupling constant, re- 
sponsible for spin flip scatterings. (The ratio of the two 
coupling constants is typically in the order of 10~ 2 .) It 
is important to note that the Hartree approximation is 
supplying a continuous BEC in the case of scalar Bose 
gases [1] , and it is related to the leading order approxima- 
tion of the 1/N expansion of the O(N) symmetric model 
in field theory [22,23]. Because of the above ambiguities a 
renormalization group approach is worked out to study the 
order of the phase transition of the homogeneous, spin-1 
Bose gas. The formulation is based on the assumption that 
the phase transition is continuous. In this case the univer- 
sal quantities (to leading order in the coupling constants) 
and the infrared (IR) behavior of the system can be calcu- 
lated from a classical field theory obtained by restricting 
the quantum fields to the zero Matsubara frequency sec- 
tor [24,8]. The main conclusion of the paper is that this 
assumption leads to contradiction which indirectly proofs 
that the transition is of first order. 
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The paper is organized as follows. In Section 2 the 
model of the homogeneous spin-1 Bose gas is mapped to 
the corresponding classical field theory by neglecting the 
nonzero Matsubara frequency components of the quantum 
fields. In Section 3 the renormalization program with min- 
imal subtraction and dimensional regularization is worked 
out up to the order of two loops for the classical field the- 
ory. In section 4 the beta function of the model is given 
up to the third power of the coupling constants and the 
renormalized free energy is derived up to one loop level. 
The critical properties of the system are discussed also in 
this section. Some summary is left to Section 5. 



2 Classical field theory 



The effective Hamiltonian of the low temperature, homo- 
geneous, spin-1 Bose gas can be written as 



V 2 are given by: 



n 

+ 



H 2 

— V# r t(aOW r (aO - ^Kx)1> r {x) 



-I* 

l - J d 3 x d 3 x% r , s s ,(x - x')*}(x)*t,(x')V s ,(x')* s (x), 



with M the mass of an atom and \x the chemical poten- 
tial. The bosonic field operator <Pjt (x) creates an atom at 
position x with spin projection r G {+1,0, —1}, and the 
operator destroys it. Automatic summation over re- 

peated spin indices is implicitly assumed throughout the 
paper. The two-particle interaction is modeled by s-wave 
scattering, with the interaction potential [15,16] 



V;, s s ,{x-x') = o (3) (x-; 



4ttH 2 
1 M 



/ \rs / \ rs 

\ / r .r g i \ / r i g r 



(2) 

where ao and a 2 are the s-wave scattering lengths in the 
total hypcrfine channel F — and F = 2, respectively. 
The matrices Vq and V 2 are the projection operators in 
the 9-dimensional tensor product space of the spin vari- 
ables projecting to the subspaces of total hyperfine spin 
and 2, respectively. The matrix Vi, projecting to the total 
hyperfine spin-1 subspace is omitted from Eq. (2), since 
it is antisymmetric in its indices and therefore does not 
appear in the Hamiltonian (1). The projection matrices 
can be obtained from the following linear equations: 



Pq + V 1 +V 2 = 1, 

-2p -r 1 + r 2 = s 1 -s 2 , 
m + v 1 + r 2 = (Si -s 2 ) 2 , 



(3a) 
(3b) 
(3c) 



where Si (i = 1,2) is the spin operator of the i.th atom. 
Using the usual spin-1 operators in the basis of S z eigen- 
vectors one can get from Eqs. (3) the projection matrices 
Vo, V\ and V 2 . Only expressing the needed two, Vo and 
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and S = 5 — 3s — s'. 



(4b) 



with R = 5 — 3r 

As shown by Ho [15] and Ohmi and Machida [16], if 
a 2 > ao the low temperature, equilibrium phase has a 
macroscopic wave function with zero net magnetization, 
called as the polar phase in analogy to the 3 He case, while 
if ao > a 2 the low temperature phase has a wave function 
with macroscopic magnetization, called as the ferromag- 
netic phase. 

In this paper we focus on the determination of the or- 
der of the phase transition. Supposing first that the possi- 
ble transitions to the polar or to the ferromagnetic phases 
of the spin-1 Bose gas are continuous phase transitions and 
restricting ourselves to universal quantities of the system, 
the problem can be mapped to a classical field theory with 
the following bare action [24,8]: 



SWM 



I 



= / d d x 



\d^;{x)d fl ip r {x) + ^m 2 ip*(x)ip r (x) 
+ Crs ' r ' s ' <p*(x)<p*,(x)<p s <(x)<p s (x) , (5) 



with tp r (x) the 3-component, complex, classical field at 
the d-dimensional position x, and is the d-dimensional 
gradient. The dimension of the system is generalized from 
3 to d for later purposes. Further on we set h = ks = 1. 
The bare mass of the field theory is denoted by m 2 = 
-2M/i, and the tensor structure of the interaction term 
takes the form of 



167r 2 a 



, \ - 167r 2 a 2 / „ \ rs 



with Ab = y/2ir/MT being the de-Broglie wavelength. 

The parameters m 2 and c rs y s > of the bare action (5) 
are easily obtained from the finite temperature action [2] 
corresponding to the Hamiltonian (1) by simply replacing 
the imaginary time dependent fields with time indepen- 
dent ones and by integrating out the imaginary time. This 
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procedure is clearly a restriction to the zero Matsubara 
frequency sector of the full field theory. The effects of the 
higher Matsubara frequency components would provide a 
physical cutoff to the field theory and change the value of 
the above parameters. These questions are not discussed 
in this paper since the order of the phase transition is not 
sensitive for the actual value of the above quantities. 

It is more convenient to introduce the matrices X = 
Vo + V2 and Y = V2 — 27^0 and to express the bare inter- 
action (6) with the help of them: 



with 



16tt 2 a + 2a 2 

16-7T 2 d2 — & 



(7) 

(8a) 
(8b) 



Note that in the case 02 > a , i.e. in the polar case c s > 0, 
while in the ferromagnetic case (a > 012) c s < 0. 



3 Renormalization up to order of two loops 

In the following we suppose that on the critical surface, 
all of the renormalized masses are zero. This assumption 
is valid in the absence of an external magnetic field and 
when the system is in a paramagnetic and non-magnetized 
phase above the (spinor) Bosc-Einstein condensation. The 
standard renormalization group program of the masslcss 
theory with dimensional regularization and minimal sub- 
traction [22,23] will be carried out for the classical field 
theory, described by the bare action (5) and interaction 
(7). The free propagator corresponding to the quadratic 
part of the bare action then reads as: 



G 



(0)rs 



(p) = S rs G {0) (p) = 6, 



p 



(9) 



which is just twice the value of the free propagator in a 
theory with real fields. 

The n-point vertex function is denoted by r^ n \pi,c), 
with c the bare interaction (7), Pi (i G {1 . . . n}) are the 
wave-numbers of the vertex function, has n spin in- 
dices. The spin indices will be omitted most of the time 
for notational simplicity. The renormalization conditions 
for the vertex functions then read as: 



r£\pi,K,g) 



Z^(g)r^( Pi ,c), (10a) 



where Z(g) is the field renormalization constant, k is the 
momentum scale, where the renormalization is made, and 
g stands for the dimensionless renormalized coupling con- 
stant, having the same index structure as the bare one 
(7). The connection between the bare and renormalized 
coupling constant is established by 



Ci]M = K e g 0ljM = K t G tjik i{g), 



(10b) 




Fig. 1. The second order divergent graph contributing to the 
2 point function. 



with e = 4 — d, go is the dimensionless, bare coupling con- 
stant and Gij t ki (g) the coupling constant renormalization 
function. The renormalization constants are expanded as 
power series of the renormalized coupling constants g%j,ki, 
according to 

Z(g) = l + b^(g,g), (11a) 
Gi 3 .ki(g) = g l3 ,u + a ( ^ kl {g,g) + a\f M (g,g,g). (lib) 

Here we have introduced the following symbolic notation 
for scalar and tensorial quadratic quantities: 



K9, 9) 



iiiji,k 1 l 1 



«2j2,K2l2 



k- 2 u9i\ji.k\l\9i'2j'2,k< 2 l- 2 j 



(12a) 

a ijki(g>g) = [ a ijki). . , , giiji,kihgi2h,k 2 h, ( 12 t>) 

respectively, and similarly for the cubic a^ kl (g , g , g) . The 
2-point and 4-point vertex functions are expanded as: 

4 2) (P. go) - P 2 % l 1 - s(2) (.90, 5o) + O(g 3 )} , (13a) 



r i^ki(Pi,go) = « £ 



90ij,M 

7(2) 



d«,fcj(0o.0b) 



+ dfAi(9o,go,go)+0(g 4 ) . (13b) 



In Eq. (13a) p 2 ■ E( 2 \g,g) is the two loop contribution to 
the self-energy. (The first order term is momentum inde- 
pendent and cancelled, since we assume that the renormal- 
ized masses are zero.) The constants and d^ are the 
corresponding 1 and 2 loop contributions to the four-point 
function. 

The divergent graph (up to the order of two loops) 
contributing to the two point function is plotted in Fig. 1. 
With the requirement that the renormalized vertex func- 
tion (10a) for n = 2 is finite, one obtains, that 



i(2)/ \ _ 9im,kl9lk,mi 

(9,9) - " g j 



(14) 



with N d = 2/(47r) d / 2 r(d/2) the usual factor appearing af- 
ter each momentum integration. r(s) is the Euler gamma 
function with argument s. 






Fig. 2. The divergent graphs contributing to the 4 pont func- 
tion at 1 loop order 
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Fig. 3. The divergent graphs contributing to the 4 pont function at 2 loop order 



The divergent graphs contributing to the four point 
function at 1 loop level are drawn in Fig. 2, while at 2 loop 
order are drawn in Fig. 3. With the requirement that the 
renormalizcd four point function [Eq. (10a) with n — 4] is 
hnite, the renormalization constants can be obtained. The 
one loop contribution is 



aW(g,g) = N d - 



» 



while the two loop contribution reads as: 
4 



'^'(9,9,9) 



e 2 l 



(2) 

(fe) 4- gV) 
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9 ■ 2(9 '/ 



e 
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(n) 
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9 {o) +9 



+ 5 W 
8(5 (e) 



9 - ■ g- + 9 (g) + 9 {h) + 2 (9 (t) + 9 U) 



(15a) 



gU) 

-9 {f) ) 

4- „(*0 



+ g {l) + 2[g^ +gW + g^ + g (p) ])] }, (15b) 



with the notations 






.(«) 
9ij,ki 


= 9im,kngmj,nh 


(16a) 


o (b) 
9ij,kl 


= gij,mngnm,kli 


(16b) 


(c) 

9ij,ki 


— gkj,mngnm,ili 


(16c) 


id) 

9ij,kl 


— gim, kngmo.npgoj.pl j 


(16d) 


(e) 

9ij,ki 


— gij,mngnm,opgpo,kl7 


(16e) 


(/) 


— gkj.mngnm.opgpo,ih 


(16f) 
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— gim,kngmj,opgpo,nl7 


(16g) 


(fc) 

yij,ki 


— gim,opgkn,pogmj,nl7 


(16h) 


a (l) 
9ij,kl 


— gij,nmgko,mpgon,pli 


(16i) 


U) 

9ij,kl 


— gkj .nmgio,mpgon.pli 


(16j) 


(fe) 

9ij,ki 


= gio,mpgoj,pngkl,nmj 


(16k) 


9ij,ki 


= gko,mpgoj,pngil,nmj 


(161) 


9ij,ki 


= gij : nmgkn,pogmp.oh 


(16m) 



ffij.fci 9kj,nm9in,po9mp,oU 

(16n) 

9ij,fci = 9io,pm9oj,np9mn,kl' 

(160) 

9r),fcZ = 9ko,pm9oj,np9mn,il- (16p) 

4 Critical properties 

The critical properties of the massless theory can be stud- 
ied with the help of the /? function, defined as 



dG 



ij,kl 



dg 



ftn 



= -eG 



ij,kl • 



(17) 



mn,op 



The /3 function can be easily expressed by inverting the 
matrix dG/dg perturbatively in g, e.g. with the help of 
iteration. The result reads as 



= -eg + 2N d 



1 ^+2{g^+g^) 



-ANi 



20 



+gW) + fl W +g& +gW +gW+2(g™ +g™ + g^+g ip) ) 



The tensorial j3 function (18) splits into two functions, 
according to 



Pij,ki = Pn{g n ,g s )Xij tk i + /3 s (g n ,g s )Y i j ! ki, 



(19) 



with X and Y defined above Eq. (7). The corresponding 
functions are: 

Pn(9n, 9s) = ~eg n + 2N d (7g 2 n + Ag n g s + Ag 2 ) 
- AN 2 (24^ + 22g 2 ngs + 3Qgng 2 + 2Qg 3 s) + 0{g 4^ (2Qa) 

and 

Ps(g n ,g s ) = -tg s + 2N d (6g n g s + 3g 2 ) 

- AN 2 (28g 2 n g s + 28g n g 2 s + g 3 s ) + 0(g 4 ). (20b) 
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Fig. 4. The flow diagram and fixed point structure at d 

(6=0). 



These functions are responsible for the flow of the renor- 
malized coupling constants under a change of scale: 



dg n ,s 



= Pn,s{9n,9s)- 



(21) 



with p being the scale, and g n ,s{p = 1) = 9n s- The IR 
behavior (p — > 0) of the system can be studied with the 
help of the IR fixed points of the (3 function (20). 

At d = 4 four real fixed points exist. The Gaussian 
one: Q, with (g n -,9s) — (0,0), and three nonphysical fixed 
points, which are absent at the one loop calculation. The 
fixed point structure and the coupling constant flow is 
shown in Fig. 4. At d = 4 — e, with < e <C 1 a new real 
fixed point, £>, emerges, which is of 0(e). The fixed point 
up to 0(e 2 ) is {g n ,9s) = (<M,0), with 



1 



<Jn 



24 



1 

-el — | e 

2N d \7 343 



O 



(22) 



This fixed point is the U(3) symmetric fixed point of Bose 
condensation. It is stable from the direction of the Gaus- 
sian one, however it is repulsive through the direction 
of g s . The fixed point structure and the flow diagram is 
shown in Fig. 5. Since there is no attractive fixed point (for 
g s 7^ 0), the trajectories "run away", which is an indica- 
tion that both the polar (when g s > 0) and the ferromag- 
netic (when g s < 0) Bose condensation is of first order. 
Such fluctuation induced first order transitions are not 
rare, e.g. the case of a real (j) 4 theory with cubic anisotropy 
[18,19]. 

The runaway of the trajectories requires an analysis 
of the free energy (thermodynamic potential) of the clas- 
sical system, or as in the terminology of field theory the 
effective action, which is the generating functional of the 
vertex functions. The free energy is obtained with the help 




-0.16 



Fig. 5. The flow diagram and fixed point structure at d < 4 
(for e = 0.4). The red lines indicate the boundary of the sta- 
bility wedge of the zeroth order contribution of the free energy 
(see below). 



of the method of steepest descent up to one loop order [22, 
23] . The zeroth order (tree-graph) contribution is the bare 
action (5): 



1 



: r n \v] = ^m 2 ip r <p r + -Cr^r's'VrVsVr'Vs', (23) 



with if being real and homogeneous. The free energy (23) 
describes a continuous phase transition at m 2 = 0. V = L d 
is the volume of the system. For m 2 < and c s > (polar 
case) the homogeneous field minimizing the potential (23) 
can be chosen as <p r = ip x (0, 1, 0) r , which has zero mag- 
netization, while for c s < (ferromagnetic case) the mini- 
mizing solution can be taken as (p r — ip x (1, 0, 0) r , which 
has maximal magnetization. With these solutions the free 
energy at tree-graph level takes the following forms: 



1 

V 
1 

V" 



polr 



„ [<p] 

<p] 



pfcrr 
I n 



™2 2 

m ip 



™2 2 

m ip 



(C„ + C s )(y2 4 . 



(24a) 
(24b) 



In the polar case (c s > 0) the tree graph thermodynamic 
potential (24a) is bounded from below only if c„ > 0, while 
in the ferromagnetic case (c s < 0) Eq. (24b) is confining 
only if c„ + c s > 0. The line of stability is therefore c„ = 
for c s > and c„ + c s = for c s < 0. As can be seen 
in Fig. 5 all interesting trajectories with g s ^ reach the 
stability boundaries of the tree graph free energy and go 
outside the region of stability. It is therefore mandatory to 
calculate the next order contribution. Following the stan- 
dard technique [22,25,26] (but with complex fields) the 
one-loop contribution of the free energy reads as 
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log 



X + m 2 + 2c. 



2 4 



+ 2 log 



(k 2 + m 2 + (c + c s )ip 2 ) 2 - c 2 ^ 4 



(25a) 



in the polar case, while 



V 1 1 1 2 



d 3 fc 



log 



(A: 2 + m 2 + 2(c„ + C;s )v5 2 



-(Cn+C s )V 4 +2 log 



m 



+ 2 log fc 2 + m 2 + ( C?1 



(c„ + c s )(p 2 



(25b) 



in the ferromagnetic case. The integrals appearing in Eqs. 
(25) are divergent. Carrying out the renormalization scheme 
renders them finite: 

\AM - f(vi) + ffo) + 2/(^3) + 2/(1/4) 



V 

both for the polar and ferromagnetic cases, with 



< o1 = 



< -- 

in the polar case and 
4° 11 -- 

ferr _ 

v 2 

ferr . 
^3 - 
, ferr 



t + 3g n ip , 

t + 9nf 2 , 

t+ (g n + 2g s )ip 2 , 
t + g n ip 2 

-t + 3(g n + g s )ip 2 , 
~-t+ (g n +g s )(p 2 , 
- 1+ (g n +g s )(fi 2 , 
t + (g n - g s )v 2 



(26a) 

(26b) 
(26c) 
(26d) 
(26c) 

(26f) 
(26g) 
(26h) 
(26i) 



in the ferromagnetic one. Here we introduced t the renor- 
malizcd, dimcnsionlcss temperature, and measured the field 
ip in units of k 1 ^ 6 / 2 , with k being the scale of the renor- 
malization. The function / appearing in Eq. (26a) is given 

by 



/(■'■) = logs + - ). 



(27) 



The sum of Eqs. (23) and (26a) gives the free energy of 
the classical system (5) in next to leading order. 

On the border of stability the free energy up to one- 
loop order (both for the polar and ferromagnetic cases) 
can be cast to the form (with neglecting terms not de- 
pending on ip): 

±r[p] = htip 2 + 2f(t + 2\g s \y 2 ). (28) 

The potential (28) describes systems with a first order 
phase transition for < \g s \ — 0(e). It is worth to note 
that in the Hartree approximation made for the quantum 
theory of the spin-1 Bose gas directly in 3 dimensions [27, 
28] , or in the equivalent lattice mean-field calculation also 
made directly in 3 dimensions [29] the phase transition 
was found to be of first order for small coupling constant 
values (at least for the ferromagnetic case). 



5 Summary 

In summary, we have studied the critical properties of 
spin-1 Bose gases with the assumption that the phase 
transition to the polar and to the ferromagnetic Bose- 
Einstein condensed phases is of second order. In this case 
the universal IR behavior of the quantum system (1) can 
be studied with the help of a classical field theory (5) 
obtained by restricting the fields to the zero Matsubara 
frequency sector. The machinery of the field-theoretical 
renormalization group then was applied to the classical 
field theory. The [3 function was determined up to the 
order of two loops in the massless theory. It was found 
that the only at least partially stable physical fixed point 
(for d < 4) is the U(3) symmetric one of Bose-Einstein 
condensation with g s = 0. However this fixed point was 
found to be repulsive towards the direction of g s . This 
indicates that all trajectories of systems with nonzero g s 
tend towards the border of the classical stability wedge 
with successive scale transformations. The free energy [on 
the boundary of the stability wedge (28)] of the classi- 
cal system (5) was determined in next to leading order in 
the method of steepest descent. It was found that both 
in the polar and ferromagnetic cases the free energy de- 
velops a second local minimum (besides the trivial one), 
which turns to be a global minimum at a certain tempera- 
ture above t = (before the second order phase transition 
sets in) . This shows that the classical system described by 
the bare action (5) exhibits a first order phase transition 
for < \g s \ = 0(e), at least near four dimensions. This 
contradicts our assumption that the polar and ferromag- 
netic Bose-Einstein condensations are continuous phase 
transitions (if they were, the classical description would 
yield also a continuous phase transition), and shows that 
the phase transition is indeed of first order both for the 
polar and ferromagnetic gases and that the jump in the 
thermodynamic quantities depend on the smaller coupling 
constant: g s . One can regard the presented approach as an 
indirect proof. 
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